q-deformed statistical-mechanical structure in the dynamics of the Feigenbaum 

attractor 
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We show that the two complementary parts of the dynamics associated to the Feigenbaum attrac- 
tor, inside and towards the attractor, form together a g-deformed statistical-mechanical structure. 
A time-dependent partition function produced by summing distances between neighboring positions 
of the attractor leads to a g-entropy that measures the fraction of ensemble trajectories still away 
at a given time from the attractor (and the repellor). The values of the g- indexes are given by the 
attractor's universal constants, while the thermodynamic framework is closely related to that first 
developed for multifractals. 



The Feigenbaum attractor, an icon of the historical 
developments in the theory of nonlinear dynamics [ij, 
is at present getting renewed attention This is be- 
cause it offers a convenient model system to explore fea- 
tures that might reflect those of statistical mechanical 
systems under conditions of phase space mixing and er- 
godicity breakdown. It therefore offers insights on the 
limits of validity of ordinary statistical mechanics. Very 
recently a thorough description with newly revealed fea- 
tures has been given 0] 01 of the intricate dynamics that 
takes place both inside and towards this famous mul- 
tifractal attractor. Here we show that these two dif- 
ferent types of dynamics are related to each other in a 
statistical-mechanical fashion, i.e. the dynamics at the 
attractor provides the 'microscopic configurations' in a 
partition function while the approach to the attractor is 
economically described by an entropy obtained from it. 
As we show below, this statistical-mechanical property 
conforms to a q-deformation Q of the ordinary expo- 
nential (Boltzmann) weight statistics. 

Trajectories within a (one-dimensional) multifractal 
attractor with vanishing Lyapunov exponent A, such 
as Feigenbaum's, show self-similar temporal structures, 
they preserve memory of their previous locations and do 
not have the mixing property of chaotic trajectories 
The fluctuating sensitivity to initial conditions has the 
form of infinitely many interlaced g-exponential functions 
that fold into a single one with use of a two-time scaling 
property 3 More precisely, there is a hierarchy 
of such families of interlaced g-exponentials; an intricate 
(and previously unknown) state of affairs that befits the 
rich scaling features of a multifractal attractor. Further- 
more, the entire dynamics is made of a family ofpairs of 
Mori's dynamical g-phase transitions 0] [ll H- 

On the other hand, the process of convergence of tra- 
jectories into the Feigenbaum attractor is governed by an- 
other unlimited hierarchy feature built into the preimage 
structure of the attractor and its counterpart repellor [j] . 
The overall rate of approach of trajectories towards the 
attractor (and to the repellor) is conveniently measured 



by the fraction of (fine partition) bins Wt^ still occupied 
at time ti by an ensemble of trajectories with initial po- 
sitions uniformly distributed over the entire phase space 
For the first few time steps the rate Wt, remains 
constant, Wj, ~ A, 1 < ii < to, to = 0(1) after 
which a power-law decay with log-periodic modulation 
sets in, a signature of discrete-scale invariance ll| . This 
property of Wt^ is explained in terms of a sequential for- 
mation of gaps in phase space, and its self-similar features 
are seen to originate in the mentioned ladder feature of 
the preimage structure The rate Wt-^ was originally 
presented in Ref. [qI where the power law exponent (f in 
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was estimated numerically. Above, h(x) is a periodic 
function with h{l) = 1, and A is the scaling factor be- 
tween the periods of two consecutive oscillations [l2|- 

We proceed now to demonstrate the connection be- 
tween the aforementioned dynamical properties. We re- 
call 1] the definition of the interval lengths or diame- 
ters dn,m that measure the bifurcation forks that form 
the period-doubling cascade sequence in unimodal maps, 
here represented by the logistic map f^{x) = I — /ix^, 
— l<x<l,{)<^<2. These quantities are mea- 
sured when considering the superstable periodic orbits of 
lengths 2", n — 1,2, 3,...; i.e., the 2"-cycles that contain 
the point a; = at /i„ < /ioo, where /ioo — 1.401155189... 
is the value of the control parameter /i at the period- 
doubling accumulation point [l3| . The positions of the 
limit 2°°-cycle constitute the Feigenbaum attractor. The 
dn^m in these orbits are defined (here) as the (positive) 
distances of the elements Xr, 



Lm, m = 

their nearest neighbors //t„ (x. 



0, l,2,...,2"-^-l, to 



(2) 



For large n, (in,o/'^n+i.o — Q!, where a is Feigenbaum's 
universal constant a ~ 2.5091. 
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Innermost to our arguments is the following compre- 
hensive property: Time evolution at /ioo from t = up 
to i oo traces the period-doubling cascade progression 
from /i = up to /ioo. Not only is there a close resem- 
blance between the two developments but also asymp- 
totic quantitative agreement. For example, the trajec- 
tory inside the Feigenbaum attractor with initial condi- 
tion xo = 0, the 2°°-supercycle orbit, takes positions Xt 
such that the distances between nearest neighbor pairs 
of them reproduce the diameters dn^m defined from the 
supercycle orbits with /i„ < /Kqo- See Fig. 1, where 
the absolute value of positions and logarithmic scales are 
used to illustrate the equivalence. This property has been 
central to obtain rigorous results for the fluctuating sen- 
sitivity to initial conditions Ct(xo) within the Feigenbaum 
attractor, as separations at chosen times t of pairs of tra- 
jectories originating close to xq can be obtained as diam- 
eters dn m, where n and m relate to t and a;o, respectively 
ii. ' 
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FIG. 1: Left panel: Absolute value of attractor positions for 
the logistic map ffiix) in logarithmic scale as a function of 
— ln(^^— /i). Right panel: Absolute value of trajectory po- 
sitions for ffiix) at with initial condition xo= in log- 
arithmic scale as a function of the logarithm of time t, also 
shown by the numbers close to the circles. The arrows in- 
dicate the equivalence between the diameters d„,o in the left 
panel, and position differences Dn with respect to xo= in 
the right panel. 

Further, the complex dynamical events that fix the 
decay rate Wt^ can be understood in terms of the cor- 
relation between time evolution at /ioo from t = up 
to < ^ 00 and the 'static' period-doubling cascade pro- 



gression from /i = up to /ioo- As shown recently 14 1 
each doubling of the period (obtained through the shift 
Mn ~^ Mn+i) introduces additional elements in the hi- 
erarchy of the preimage structure and in the family of 
sequentially-formed phase space gaps in the finite period 



cycles. The complexity of these added elements is simi- 
lar to that of the total period 2" system. Also, the shift 
Mn ^ A^n+i increases in one unit the number of undula- 
tions in the transitory log-periodic power law decay found 
for the corresponding rate Wn,ti of approach to the 2^' 
supercycle attractor 



14| 



As a consequence of this we 
have obtained detailed understanding of the mechanism 
by means of which the discrete scale invariance implied 
by the log-periodic property in Wt^ = lim„^oo W^n,ti 
arises. What is more, the rate Wt^ , at the values of time 
for period doubling, can be obtained quantitatively from 
the supercycle diameters ci„_m, that is [1], Wt^ — IS. Zt, 
t = ti — tQ, and 



271- 



Zt = 



dn- 



i = 2" 



n= 1,2,3, , 



(3) 



m=0 



Eq. ([31) is an explicit expression equivalent to the nu- 
merical procedure followed in Ref. [l^ by the use of 
the triadic cantor set construction of the Feigenbaum at- 
tractor to evaluate the power law exponent (p, and from 
which the value for (p = 0.800138194 is reported. 

Now, to reveal the aforesaid statistical-mechanical 
structure we identify the decay rate Zt as a partition 
function. From this perspective the diameters dn,m are 
configurational terms and we go forward to determine 
their time dependence through their indexes n and to. 
The dn rn scalc with n for to fixed as 



dn.r. 



„ — 1 2 3 



(4) 



the dim are universal constants obtained, for instance, 
from the finite jump discontinuities of Feigenbaum's tra- 



jectory scaling function a{y) = lim^ 



^00 ^n,m+l 



/dn,, 



y = lim„^oo to/2" [i[. The largest discontinuities of a{y) 
correspond to the most sparse and the most crowded 
regions of the multifractal attractor, and for these we 
have, respectively, c?„^o — a"""^"^ and dn,i — 
(The 1st diameter do,o — 1 and the equality in Eq. 
(HI) is rapidly reached as n increases). With use of the 
identity ^""+1 = (1 + /? = 2""^ - 1, the 

power law in Eq. ([4]) can be rewritten as a g-exponential 
(exp,(x) = [l-(g-l)x]^i/(9-i)),i.e., 



where Qm 

-1 



1 



dn,m =i exp^^^(-/3i/„), 

lna,„/ln2, and f3 — t - 



(5) 



I/, 



1 = 

2"^^ — 1. Likewise, the partition function Zf ~ t^*^ (or 
Zt ~ with ip = lne/ln2 and t = 2"-^, can be 

written as 



?t ~ expQ(-/3(/3), 



(6) 



Q = 1 + (^-1 and again p = t-l = 2"-^ - 1. 

Our main point in this Letter becomes apparent when 
Eqs. (O and © are used in Eq. ([3]), to yield 



expQ{-(3(p) ~ ^exp^^^(-/3j/„). 



(7) 
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Eq. ([7]) resembles a basic statistical-mechanical expres- 
sion where the quantities in it play the following roles: 
(3 an inverse temperature, a free energy (or the prod- 
uct s = —Pf a Massieu thermodynamic potential, or 
entropy), and the configurational energies. However, 
the equality involves g-deformed exponentials in place of 
ordinary exponential functions that would be recovered 
when Q = Qm = 1- It is worth noticing that there is 
a multiplicity of q-indexes associated to the configura- 
tional weights in Eq. ([7]), however their values form a 
well-defined family [6[ determined by the discontinuities 
of Feigenbaum's function a. To substantiate the useful- 
ness and appropriateness of this identification we present 
in the remaining part of this Letter: 1) A 'mean field' 
evaluation of Zt and a thermodynamic interpretation of 
the time evolution process. 2) The relationship of Zt with 
the familiar partition function developed for the descrip- 
tion of multifractal geometry. 3) A crossover to q = 1 
ordinary statistics. And 4) a commentary on the im- 
plications of our results for the o-deformed generalized 
statistical mechanical formalism [a. 




X 

FIG. 2: Sector of the bifurcation tree for the logistic map 
f^{x) that shows the formation of a Pascal triangle of diame- 
ter lengths according to the scaling approximation expalined 
in the text, a ~ 2.5091 is the pertinent universal constant. 

Akin to a mean field approximation we assume that, 
for a given value of n, e.g., n = 3, the diameters dn,m that 
are of comparable lengths have equal length and this is 
obtained from those of the shortest or longest diameters 
via a simple scale factor; e.g., ^3^3 = ^3^2 = <^~^d3,a — 
ad^^i. This introduces some degeneracy in the lengths 
that propagates across the bifurcation tree. See Fig. 2 
[l^ . Specifically, the dn,m scale now with increasing n ac- 
cording to a binomial combination of the scaling of those 
diameters that converge to the most crowded and most 
sparse regions of the multifractal attractor. This is to 



consider that the 2""-*^ diameters at the n-th supercycle 
have lengths equal to q;~''"~^~''q;~^' and occur with mul- 
tiplicities ("7^)' where ^ = 0, 1, n — 1. As seen in Fig. 
2 the diameters form a Pascal triangle across the bifurca- 
tion cascade. The partition function can be immediately 
evaluated to yield 



(8) 

t = 2"-i. We obtain e = {a~^ + a'^) = 1.7883, 
(p — 0.8386, and Q = 2.1924, a surprisingly good ap- 
proximation when compared to the numerical estimates 
(p = 0.8001, and Q = 2.2498 of the exact values. Un- 
der this approximation all the indexes qm in Eq. ([7]) are 
equal, g,„ = 5 = 1 -I- v~^, v — In a/ In 2, and Eq. ([7]) 
becomes 

ri-l 

expQ(-/3(^) = - 1' exp,(-/3i/), (9) 

where D,{n —\,l) = a~' ("7^)' Thermodynamically, the 
approach to the attractor described by Eqs. ([7]) or (l9|) is 
a 'cooling' process — > 00 in which the free energy (or 
energy) Lp is fixed and therefore the entropy s = —f3ip is 
linear in /3. It is illustrative to define an 'energy land- 
scape' for the Feigenbaum attractor as being composed 
by an infinite number of 'wells' whose equal-valued min- 
ima at /3 — > 00 coincide with the points of the attractor 
on the interval 1] [H. When f3 = 2""! - 1, n 

finite, the wells merge into 2"~^ intervals of widths equal 
to the diameters dn,m- 

As it is well-known the geometric properties of multi- 
fractals conform to a statistical mechanical framework, 
the so-called thermodynamic formalism [l3|. The parti- 
tion function devised to study their properties, such as 
the spectrum of singularities f{a) ^13] . is 

M 

Z(T,q)^^p;,C^ (10) 

m 

where the Im (in one-dimensional systems) are M dis- 
joint interval lengths that cover the multifractal set and 
the Pm are probabilities assigned to these intervals. The 
usual procedure consists of requiring that 2'(t, q) nei- 
ther vanishes nor diverges in the limit Im for all 
TO (and consequently M — > 00). In this case the expo- 
nents T and q define a function r(q) from which f{a) 
is obtained via Legendre transformation [l^. When the 
multifractal is an attractor its elements become ordered 
dynamically, and for the Feigenbaum attractor the tra- 
jectory with initial condition xq = generates sequen- 
tially the positions that form the diameters, producing 
all diameters dn,m for n fixed between times t = 2"~^ 
and t = 3 • 2"~^. Since the diameters cover the attractor 
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it is therefore natural to choose the covering lengths at 

fn) 

stage n to be /„/ = c?n,m and to assign to each of them 
the same probability p™'' — 1/2. For instance, within the 
two-scale approximation to the Feigenbaum multifractal 
[l^, = a-("-i-'^-)a-2fe, the condition Z(r,q) = 1 re- 
produces Eq. dH) when p^^^ = = 2-"-+i, with r = 1 
and q ~ —ip. It should be kept in mind that the 'static' 
partition function Z(t, q) is not meant to distinguish be- 
tween chaotic and critical (vanishing A) multifractal at- 
tractors as we do here. As we emphasize below, it is 
the functional form of the link between the probabilities 
Pm ^ and actual time t that determines the nature of the 
statistical mechanical structure of the dynamical system. 

The recursive method of backward iteration of chaotic 
maps provides a convenient way for reconstructing multi- 
fractal sets and obtaining their underlying statistical me- 
chanics [T3|- A chaotic unimodal map has a two- valued 
inverse and given a position a; = a;„ a binary tree is 
formed under backward iteration, so there are 2" initial 
conditions xq for trajectories that lead to Xn. Since the 
Lyapunov exponent is positive A > 0, for large n lengths 
expand under forward iteration according to I exp(An) 
and contract under backward iteration as I ^ exp(— An). 
We can define, as above, a set of covering lengths (5„,„i = 
exp(— Am'^), where A^ is a local Lyapunov exponent with 
-^m — > A as n — !■ cx) and where m relates to the initial con- 
dition xq. Use of these lengths in a partition function like 
that in Eq. ^ gives 

exp(-/3$) =^exp(-/3A,„), (11) 

m 

where now (3 ^ n. Recalling Pcsin's theorem [13], <& 
is clearly related, for large n, to the Kolmogorov-Sinai 
entropy. The crossover from g-deformed statistics to 
ordinary q = 1 statistics can be observed for control 
parameter values in the vicinity of the Feigenbaum at- 
tractor, > fXoo: when the attractor consists of 2" 
bands. The Lyapunov coefficient A of the chaotic at- 
tractor decreases with A/i = fi— fjLoo as A oc 2^" ~ A/j,'^, 
K, = In 2/ ln(5i?(<r), where 5f is the Feigenbaum constant 
that measures the rate of development of the bifurca- 
tion tree in control parameter space [ij. The chaotic or- 
bit consists of an interband periodic motion of period 
2" and an intraband chaotic motion. The expansion 
rate Yl\=^ ^'^ \dffi{xi)/dxi\ fluctuates with increasing am- 
plitude as In t for i < 2" but converges to a fixed number 
that grows linearly with t for t ^ 2" Q . This translates 
as dynamics with q ^ 1 for < < 2" but ordinary dynamics 
with g = 1 for t > 2" 0. 

We have shown that there is a statistical-mechanical 
property lying beneath the dynamics of an ensemble 
of trajectories en route to the Feigenbaum attractor 
(and repellor). The fraction of phase space still occu- 



pied at time t is a partition function Zt made up of 
q-exponential weighted configurations, while Zt itself is 
the g-exponential of a thermodynamic potential func- 
tion. This is a clear signature of g-deformation of or- 
dinary statistical mechanics, and, to our knowledge, it 
is the first bona fide concrete instance (anticipated or 
not in the form presented here) where arguments can 
be made explicit and rigorous. There is a close parallel 
with the thermodynamic formalism for multifractal sets, 
but it should be emphasized that the departure from the 
usual exponential statistics is dynamical in origin, and 
due to the vanishing of the (only) Lyapunov exponent. 
Our results suggest a novel variant for the theory of large 
deviations, the mathematical articulation of statistical 
mechanics [l^. 
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